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I.  INTRODUCTION 


An  important  parameter  used  In  gun  Interior  ballistic  coaputatlons  Is  the 
propellant's  lapetus,  the  energy  contained  In  the  hot  chaaber  gases  produced 
by  one  graa  of  the  propellant  when  It  Is  burned  to  coapletion  at  a  specified 
loading  density  (coaaonly  0.2  g/cc).  The  coaputatlon  of  the  lapetus  for  the 
general  case  is  straightforward  only  In  principle;  In  practice  aany 
coaplications  arise. 

The  most  obvious  one  is  essentially  nuaerical.  In  a  typical  case,  the 
chamber  gas  contains  soae  50-60  species  present,  only  a  few,  usually  less  than 
10,  have  any  appreciable  concentration.  As  teaperatures  and  pressure  change, 
the  identities  of  these  10  species  also  change.  The  concentrations  of  all  of 
the  species,  the  temperature,  and  the  pressure  are  related  to  one  another  by 
complex  non-linear  equations.  Overall,  the  laws  of  conservation  of  mass 
energy  aust  be  obeyed. 

This  problem  Is  complicated  enough  for  a  mixture  containing  only  gases; 
the  presence  of  one  or  more  condensed  phases  (e.g.,  potassiua  hydroxide, 
carbon,  alunlnua  oxide)  makes  it  even  anre  coaplicated.  The  concentrations  of 
all  of  the  gaseous  species  ,  but  not  the  condensed  phases,  obey  the  familiar 
"mass  action-  equations,  and  hence, the  latter  Bust  be  treated  separately. 

The  general  case  was  not  atteapted  until  the  developaent  of  electronic 
computers  aade  it  feasible.  The  need  for  such  coaputatlons  became  a  large 
national  enterprise;  aany  computer  codes  were  written  for  this  purpose.  Since 
rocket  aotors  usually  operate  below  35  MPa,  these  codes  all  essuaed  the  Ideal 
gas  law,  pV-nKT. 

Entirely  different  conditions  are  encountered  in  the  detonation  of 
condensed  explosives.  Here,  the  pressures  exceed  hundreds  of  kllobars.  The 
type  of  gas  law  that  aust  be  used  under  these  conditions  Is  controversial,  and 
is  outside  the  scope  of  this  paper. 

The  initial  chaaber  conditions  for  guns  typically  have  gas  densities  cf 
0.2  g/cc  or  higher;  at  chese  densities,  ideal  gas  pressures  are  10  to  20Z  too 
low,  so  a  correction  -»<st  be  aade.  A  way  of  doing  this  was  first  pointed  out 
by  Corner.*  His  nethed  could  and  was  iapleaented  using  mechanical  desk 
calculators. 

The  first  electronic  coaputer  prograa  embodying  Corner's  Ideas  was 
written  circa  1950  at  the  Ballistic  Research  Laboratories  (now  the  Ballistic 
Research  Laboratory),  but  it  received  little  notice.  A  report2  containing 
results  computed  with  this  program  oid,  however,  receive  wide  circulation. 


JJ.  Comer,  Theory  of  the  Interior  Ballietioe  of  Guns.  Wiley,  New  York 
(1950). 

2P.  G.  Baer  and  K.  R.  Bryson,  "Tables  of  Confuted  Thermodynamic  Properties 
of  Military  Sun  Propellents,"  BEL  Memorandum  Report  1338  (1961).  AD-258-644 


More  recently,  the  present  writer  took  advantage  of  the  nodular  design  of 
the  TIGER  program,  which  permits  an  alaost  unlimited  number  of  different 
gaseous  equations  of  state  to  be  Inserted  into  it  with  relative  ease,  and 
programmed  Corner's  method  for  it.  The  resulting  code,  called  BLAKE,  is  now 
used  at  several  US  military  installations,  including  AMCCOM-Dover,  AMCCOM  BRL, 
and  the  Naval  Surface  Weapons  Laboratories  at  Indian  Head  and  Dahlgren, 

Even  more  recently,  Volk  and  Bartheldt^  published  a  description  of  a  code 
that  accomplishes  the  same  task.  While  there  are  of  course  some  differences, 
it  appears  that  their  approach  was  qualitatively  similar. 

As  a  result  of  Dr.  F.  Volk's  visit  to  the  BRL  in  June  1977,  it  was 
suggested  that  a  round  robin  be  conducted  between  the  various  available 
codes.  Dr.  Volk  agreed  to  compare  the  various  FRG  codes,  and  this  writer 
undertook  to  do  the  same  for  the  US  codes, 

II.  OTHER  US  CODES 

If  one  thermodynamics  code  is  to  be  singled  out  and  labeled  the  "best , " 
surety  it  must  be  the  NASA/ Lewis  code,  CEC71,  first  written  almost  30  years 
ago  and  its  revised  and  extended  successor,  TRAN72.  It  offers  trouble- 
free  operation,  a  wide  range  of  applications,  and  the  JANNAF  library  of 
reaction  products.  The  only  drawback  to  this  marvelous  program  is  that  it  has 
the  ideal  gas  equation  of  state  Inextricably  woven  into  it. 

The  calculations  presented  for  the  NASA/Lewis  program  were  performed  by 
R.  A.  Flfer,  Ignition  and  Combustion  Branch,  Interior  Ballistics  Division, 

BRL. 

There  are,  of  course,  a  large  number  of  thermodynamics  programs 
available,  of  which  three  have  been  widely  applied  to  propellants.  These  are 
TIGER,  PEPO,  AND  NOTS.  There  was  nothing  to  be  learned  by  using  the  original 
TIGER  program  for  the  present  exercise,  since  this  program  is  Incorporated 
into  BLAKE,  as  already  noted.  The  revised  version  of  TIGER  was  not  available 
for  this  project. 

PEPO  (an  acronym  derived  from  Propellant  Evaluation  Program)  has  been 
extensively  used  in  many  US  Navy  laboratories.  It  is  limited  to  ideal  gases 
and  cannot  do  constant  volume  calculations  and  so  was  not  further  considered. 


•V.  Volk  and  H.  Bartheldt.  Propellants  and  Explosives.  Vol  1,  7-  14,  (1976). 

4 (a)  S.  Gordon  and.  F.  Zeleznik,  "A  General  IBM  704  or  7090  Computer  Program 
for  Computation  of  Chemical  Equilibrium  Compositions. ..."  NASA  TN  D-14S4 
(October  1962). 


(b)  R.  Soehla  and  B.  McBride,  " Fortran  IV  Computer  Program  for  Calculation  of 
Thermodynamic  and  Transport  Properties  of  Complex  Chemical  Systems, "  NASA  TN 
D-70S6.  ( January  1973). 


NOTS  (which  was  first  written  at  the  US  Naval  Ordnance  Test  Station, 
hence  Its  name)  Is  a  descendant  of  PEPO.  It,  too.  Is  Halted  to  the  ideal  gas 
state,  but  it  can  do  constant-voluae  calculations.  NOTS  Is  extensively  used 
at  the  Naval  Weapons  Center  at  China  Lake,  California.  Previous  lnforaal 
comparisons  between  it  and  BLAKE  indicated  that  the  two  codes  agreed  well  for 
the  ideal-gas  case.  Efforts  at  including  NOTS  in  the  present  study  were 
unsuccessful  owing  to  a  high  priority  project  at  China  Lake.  The  BLAKE  code 
has  been  described  elsewhere. 

III.  INPUT 

The  well-known  Encyclopedia  of  Propellants^3  was  used  as  the  source  for 
composition  data  on  five  different  US  military  powders.  Their  compositions 
were  normalized  so  as  to  add  to  100  percent,  including  residual  solvent  and 
additives,  if  any.  (Bear  in  mind  that  US  military  powder  description  sheets 
do  not  usually  present  component  information  directly  as  percentages  of  the 
final  mix.)  The  values  chosen  are  given  in  Table  i. 

These  compositions  were  chosen  to  satisfy  several  requirements. 

Obviously,  they  should  be  typical  of  modern  gun  propellants.  Additionally,  it 
is  desirable  that  they  cover  a  range  of  temperatures  and  energies.  For  double 
base  propellants  of  the  1)1)68  used,  energies  and  flame  temperatures  are 
usually  closely  correlated.  If  further  testing  of  the  codes  is  required,  a 
nltramine  composition  for  which  this  correlation  does  not  hold  should  be 
included. 

Compared  to  either  NASA/Lewis  or  the  NOTS  codes,  BLAKE  runs  more  slowly 
with  compositions  for  which  solid  products  are  formed  at  chamber  conditions. 
Also,  convergence  problems  arise  in  such  cases.  Composition  M15  was  picked  to 
represent  this  class  of  propellants.  The  BLAKE  calculation  for  it  converged 
readily  as  soon  as  the  list  of  possible  gaseous  products  was  trimmed  to  29, 
the  maximum  permitted  at  any  one  time. 

Dr.  Volk  and  this  writer  had  earlier  agreed  that  the  BRL  would  furnish 
all  of  the  heats  of  formation  for  the  initial  ingredients.  These  values  are 
given  in  Table  2.  Both  groups  rely  strongly  on  the  various  standard 
compilations  of  thermochemical  data,  especially  the  books  by  Cox  and 
Pilcher,”1’  and  by  Stull,  Westrum,  and  Sinke,  so  it  is  unlikely  that  there 
would  have  been  serious  disagreements  on  any  ingredient  other  than 


i,E.  Freedman,  "BLAKE  -  A  Thermodynamics  Code  Based  on  TIGER:  Users '  Guide 
and  Manual,"  ARBRL-TR-2411  (July  1982).  ADA-121-259. 

60 Basil  T.  Federoff  and  Oliver  E.  Sheffield ,  Encyclopedia  of  Explosives  and 
Related  Items.  PATR  2700,  Vol  2,  pp  t-34,35  (1962). 

e^J.  V.  Cox  and  G.  Pilcher,  Thermochemistry  of  Organic  and  Oraano-metdllic 
Compounds.  Academic  Press,  London  (1970). 

7D.  R.  Stull,  E.  F.  Westrum,  Jr.,  and  G.  C.  Sinke,  The  Chemical 
Thermodynamics  of  Organic  Compounds,  Wiley,  Rea  York  (1969). 


TABLE  1.  SOME  TYPICAL  US  MILITARY  POWDERS 


TABLE  2.  NAMES.  FORMULAS.  AND  KEATS  OF  FORMATION 


and  molecular  weight  =  162.143  +  44.998 


nitrocellulose.  For  this  substance,  the  BRL  relies  on  the  values  published  by 
Jessup  and  Prosen, ®  while  ICT  relies  on  Medard.** 

IV.  THE  WEST  GERMAN  ROUND  ROBIN 

Details  of  the  West  Genian  Round  Robin  will  be  handled  by  Dr.  Volk, 1(1 
who  nay  publish  the  results  separately. 

V.  THE  BLAKE  NASA/LEWIS  COMPARISON 

This  comparison  was  nade  on  the  basis  of  ideal  gas  behavior  only,  for 
the  reason  noted  above.  The  results  are  shown  in  Table  3.  They  are  obviously 
in  excellent  agreement.  An  additional  indication  of  the  close  agreement 
occurs  with  M15,  where  both  codes  show  the  formation  of  liquid  AljO*.  BLAKE 
gives  the  nole  fraction  of  this  product  as  0.000124;  NAS  A/ Lewis ,  0.0001236. 

The  only  significant  difference  between  the  two  sets  of  computations 
occurs  with  the  ratio  of  heat  capaclt.'es,  ganaa,  and  the  heat  capacity  at 
constant  pressure,  Cp.  This  will  be  discussed  in  detail  below. 

VT.  COMPARISON  BETWEEN  BLAKE  AND  THE  ICT  CODE 

Unfortunately,  a  comparison  between  the  two  codes  for  the  ideal  gas  ease •' 
cannot  be  made,  since  Dr.  Volk*1  sent  results  only  for  the  real  gas  case. 

The  results  for  the  two  codes  are  shown  in  Table  4. 

The  agreement  between  the  calculated  temperatures  is  such  that  it  appears 
probable  that  the  two  codes  would  agree  well  for  the  ideal  gas  case;  still,  if 
would  be  desirable  to  make  the  comparison  under  those  conditions. 

Overall,  the  agreement  is  good.  At  higher  loading  densities,  the  BLAKE 
results  are  lower  than  the  ICT  results  (or  the  latter  are  higher).  The  former 
case  is  the  more  likely,  since  BLAKE  computes  the  third  virial  coefficient  on 
the  assumption  of  hard  spheres  with  a  radius  0.81  times  the  Lennard-Jones 
radius,  as  suggested  by  Hirschfelder.  The  ICT  code  calculated  both  second 
and  third  virial  coefficients  from  either  the  Stockmayer  or  the  Lennard-Jones 
potentials,  depending  on  whether  the  molecule  is  polar  or  not. 


eR.  S.  Jessup  and  E.  J.  Proeen,  J.  Research  NBS  44  ,  387  (1950). 

9L.  Hedard,  Hem.  Artillerie  Francaiee  (1954). 

I0F.  Volk,  private  communication  to  E.  Freedman,  dated  22  December  1977. 

Volk,  private  communication  to  E.  Freedman,  dated  19  January  1978. 

12J.  0.  Hirechfelder,  C.  F.  Curtise,  and  R.  B.  Bird,  Holecular  Theory  of 
Fluids ,  John  Wiley,  Few  Fork,  p  157  (1954). 
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A  mixture  of  reacting  gases  has  two  heat  capacities.  The  first  of  these 
is  merely  the  mole-fraction-weighted  sum  of  the  individual  heat  capacities  of 
the  components  of  the  mixture.  The  other  heat  capacity  consists  of  the  frozen 
value  plus  an  additional  term  that  arises  from  the  contribution  of  the 
shifting  equilibrium.  This  additional  term  is  roughly  proportional  to  the 
square  of  the  enthalpy  of  reaction.  The  matter  has  been  discussed  by  a  number 
of  authors;  see,  for  example,  Hirshcfelder  et  ml.  For  an  application  to 
propellants,  see  Fifer  and  Cohen.13 

The  "transport*  option  of  the  modified  NASA/ Lewis  program4*5  computes  both 
heat  capacities.  The  results  for  a  cool  propellant  (K1A1)  and  a  hot  one  (M9) 
are  shown  in  Table  5. 

The  agreement  among  the  various  heat  capacities  and  hence  the  gammas  for 
M1A1  is  good  for  all  of  the  codes.  With  H9,  however,  only  the  frozen  heat 
capacities  agree;  the  equilibrium  value  is  about  50Z  larger  than  the  frozen 
one. 

The  heat  capacity  itself  is  not  an  extremely  important  quantity  in 
interior  ballistics,  but  gamma  is.  As  C  increases,  Cy  also  Increases,  so 
gamma  is  not  a  trivial  natter.  Apparently, there  are  enough  undetermined 
factors  remaining  in  the  treatments  of  Interior  ballistics  that  this  effect  is 
not  yet  significant.  Nevertheless,  it  deserves  explicit  treatment. 

VII.  CONCLUSIONS 

THE  NASA/ Lewis  code  and  the  BLAKE  code  agree  virtually  perfectly  on  the 
calculation  of  temperature,  pressure,  and  related  quantities  (except  heat 
capacity  and  gamma)  for  hot  propellant  gases  considered  to  be  ideal  gases. 

The  BLAKE  code  and  the  ICT  code  are  In  very  good,  but  not  perfect 
agreement. 

The  area  of  disagreement  between  BLAKE  and  the  ICT  code  is  very  likely 
the  differences  in  the  correction  for  non-ideal  gas  effects. 

It  Is  not  obvious  yet  that  these  differences  are  great  enough  to  require 
further  work. 

The  use  of  frozen  heat  capacities  (and  hence  gammas)  in  gun  codes  should 
be  closely  examined. 


A.  Fifer  and  A.  Cohen,  "High  Temperature  Thermal  Conductivity  and  Beat 
Transfer  -in  the  Gun  Propellant  and  Primer  Combustion  Gates,”  Proceedings  of 
the  Hth  MHBAF  Cmbustion  Nesting  CBIA  Publication  292,  December  1977, 
Volume  II,  p  363, 
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7..  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Organization 


City,  State,  Zip 


(Remove  this  sheet  along  the  perforation,  fold  as  indicated,  staple  or  tape 
closed,  and  nail.) 
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Director 

OS  Amy  Ballistic  Research  Laboratory 
ATTO:  DRXBR-OD-ST 

Aberdeen  Proving  Ground,  HD  21005-9989 


BUSINESS  REPLY  MAIL 

FUST  CLASS  PERMIT  NO  12062  WfcSHNGT0H,PC 
POSTAGE  WILL  BE  WUO  BY  DEPARTMENT  OF  THE  AJ*4Y 


